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ABSTRACT 

We find, from our study of binary spiral galaxies in the Sloan Digital Sky Sur- 
vey Data Release 6, that the relative orientation of disks in binary spiral galaxies 
is consistent with their being drawn from a random distribution of orientations. 
P4 For 747 isolated pairs of luminous disk galaxies, the distribution of 0, the angle 

between the major axes of the galaxy images, is consistent with a uniform distri- 
bution on the interval [0°,90°]. With the assumption that the disk galaxies are 
oblate spheroids, we can compute cos /3, where /3 is the angle between the rotation 
axes of the disks. In the case that one galaxy in the binary is face-on or edge-on, 
the tilt ambiguity is resolved, and cos (3 can be computed unambiguously. For 94 
isolated pairs with at least one face-on member, and for 171 isolated pairs with 
at least one edge-on member, the distribution of cos (3 is statistically consistent 
with the distribution of cosi for isolated disk galaxies. This result is consistent 
with random orientations of the disks within pairs. 

Subject headings: galaxies: fundamental parameters — galaxies: photometry — 
galaxies: statistics 



INTRODUCTION 



The structure of a spiral galaxy is closely tied to its angular momentum distribution. 
Thus, understanding galaxy formation requires knowledge of how galaxies acquire their an- 
gular momentum. The tidal torque model for the acquisition of angular momentum by 
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galaxies (lHoyldll949t lPeebleslll969l ) states that protogalaxies acquire their angular momen- 
tum through the tidal interaction between the protogalaxy and the surrounding nonuniform 
matter distribution. Spiral galaxies are particularly useful for testing how galaxies acquire 
angular momentum. Spiral galaxies contain relatively cold, thin disks of stars, gas, and 
dust, supported primarily by rotation rather than velocity dispersion. These thin disks are 
unlikely to have suffered a major merger (or large numbers of minor mergers) since their for- 
mation; thus, the orientation of a disk's spin angular momentum is unlikely to have changed 
drastically since the material that formed it was initially torqued. 

Binary spiral galaxies are particularly useful test cases for the acquisition of angular 
momentum. A "binary spiral galaxy" can be defined as a pair of spiral galaxies that are 
relatively close to each other, but are much farther away from other galaxies of comparable 
or greater mass. Our galaxy and M31, using this general definition, can be thought of as 
a binary spiral galaxy. They are separated by ~ 0.8 Mpc, and M33, the next brightest 
galaxy in the Local Group, is ~ 2 magnitudes fainter than our galaxy or M31. If a pair of 
spiral galaxies is a totally isolated system, with zero net angular momentum initially, then 
conservation of angular momentum states that the net angular momentum remains zero after 
the two galaxies have spun each other up. If the net angular momentum of the binary system 
is, in fact, zero, then the relative orientation of the internal angular momentum vectors of 
the two galaxies can range from anti-parallel, if the two galaxies are on radial orbits relative 
to each other, to parallel, if the orbital angular momentum is large. 



In the real universe, binary spiral galaxies are never completely isolated. iGott &: Thuan 



( 119781 ). for instance, in their study of the angular momentum of the Local Group, conclude 
that most of the angular momentum of our galaxy was produced by M31, but that the 
tidal torquing by nearby groups, such as the Maffei group {d ~ 3 Mpc) and M81 group 
{d 3.6 Mpc), was not negligible. Thus, we might expect the distribution of the angle a 
between the internal angular momentum vectors of the two galaxies in a pair to depend 
on the number of massive galaxies within several megaparsecs of the pair. In addition, 
the distribution of a will also depend on the evolutionary history of binary spiral galaxies 
since they initially attained their angular momentum by tidal torquing. For instance, the 
merge r timescale for a bound pair of galax ies depends on the ellipticity e of their mutual 



orbit JLacev fc Colelll993l : [jia^g et al.ll2008h . with more nearly radial orbits having a shorter 



merging timescale. Thus, binary spiral galaxies with radial orbits will more rapidly merge 
to form an elliptical galaxy. In the case that the binary spiral galaxies are initially isolated, 
this will result in preferential destruction of the binaries whose internal angular momenta 
are anti-parallel. 



Computer simulations have provided additional insight into the alignment of galaxy 
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spins. In n-body A cold dark matter simulations, for instance, massive dark halos (M > 
3 X IO^^Mq) have spins that are p referentially oriented in a direction perpendicular to the 
mass distribution ( jPaz et al.l |2008[ ). and have proje cted major axes that are aligned with 



the surrounding galaxy distribution (iPaz et al.l 120 111 ): these alignments persist on scales up 
to ~ 30 Mpc. It should be kept in mi nd, though, that a dark halo is not necessarily a 



simple structure; ISchneider. et al.l ( l201ll ) show that in the Millennium-2 n-bod y simulation 



25% of all halos have nearly perpendicular major axes at large and small radii. iHahn et al. 



(120101 ). in a hydrodynamic adaptive mesh refinement simulation that permits disks to form 
by dissipation, find that at 2; = 0, the spin of stellar and gaseous disks are well aligned with 
the spin of the inner dark halo; however, the alignment between disk spin and the spin of 
the entire host halo is significantly weaker. 

Many studies have been done of the relative orientation of the spin vectors of neighboring 
spiral galaxies. Three different methods, using increasing amounts of information for each 
binary spiral galaxy, may be distinguished. 

The first method for studying the relative orientation of disk galaxies uses only the 
position angle of the apparent major axis of each galaxy. If a galaxy is approximated as 
a rotationally flattened oblate spheroid, then its apparent major axis is at right angles to 
the projection of its spin axis onto the plane of the sky. The angle between the apparent 
major axes of the two galaxies in a binary spiral galaxy, defined so that 0° < (p < 90°, is 
thus the angle between the projections of the spin axes. If the spin axes of the two spiral 
galaxies tend to be parallel (implying either parallel or anti-parallel spin vectors), then the 
distribution of will be weighted toward = 0°; if the spin a xes tend to be perpendicular, 
then the distribution of will be weighted toward = 90° (ISofuelll992l). Employing this 



metho d, ISharp et al.l ( 1l979l ). with a sample of 57 pairs of spirals, and ICervantes-Sodi et al. 
(l2010l ). with 218 pairs of spirals, found distributions for that were statistically consistent 



with being uniform from 0° to 90°; th ese results ar e consistent with, but do not demand, 
uncorrected spin vectors. By contrast, ISofud (119921 ) . using a sample of 390 pairs of spirals, 
found a significant excess of pairs with ~ 90°, indicating a tendency for spin axes to be 
perpendicular. 

The second method for studying the relative orientation of spiral galaxies uses the 
apparent axis ratio q of each galaxy image in addition to the angle between their apparent 
major axes. If each galaxy is approximated as an oblate spheroid of intrinsic short-to-long 
axis ratio 7, then the inclination i of a galaxy with apparent axis ratio q is given by the 
usual relation 



cos? 



q^ 



'Jl 

7^ 



1/2 



(1) 
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Now consider a pair of disk galaxies, with inclinations ii and Z2, that have an angle (f) between 
their apparent major axes. If the angular separation of the galaxies is small, then the angle 
P between their spin axes, defined so that < /3 < 90°, is given by the relation 



cos /3 = I sin 2i sin ^2 cos ± cos ii cos ^2 1 



(2) 



The two-fold ambiguity in cos /3 is due to the tilt ambiguity of the individual galaxies. (If 
the apparent major axis of a disk galaxy lies in the east-west direction, for example, we 
don't generally kno w whether the northern or the southern half of the disk is closer to us.) 
Using this method, iFlinl ( 119931 ). with a sample of 586 pairs of galaxies, found that their spin 
axes tend to be parallel. However, this method is intrinsically unable to distinguish between 
parallel or anti-parallel spin vectors. 

The third method for studying relative orientation of spiral galaxies uses additional 
information to resolve the tilt ambiguity and spin ambiguity of disk galaxies. The tilt 
ambiguity ("Is the northern or southern half of the disk closer to us?") can be resolved 
by looking at high-resolution images of a disk galaxy and assuming that dust lanes lie 
on the outer edge of spiral arms. The spin ambiguity ("Is the disk rotating clockwise or 
counterclockwise from our point of view?" ) can be resolved, if the galaxy is close to face-on, 
by looking at the chirality of the disk's spiral arms and assuming that spiral arms are trailing. 
Either of these two pieces of information - location of dust lanes and chirality of spiral arms 
- can be replaced by spectroscopic information about which half of the disk is redshifted 
relative to the galaxy's nucleus and which is blueshifted. The additional information lets 
us find the angle a betwee n the s pin vectors of the galaxies, defined so that < a < 180°. 
Using this method, iHeloul ( 1l984j ). with a sa mple of 31 pairs of galaxies, found that spin 



vectors tend to be anti-parallel 



a 



180°). lOosterlod (Il993f), using asample of 40 pairs 



Pestana fc Cabreral (l2004r) used a sample 



found that their spin vectors were uncorrelated. 

of 46 pairs of galaxies, largely drawn from the study of Oosterlool ( 19931 ). but discarded the 
information that resolved the spin ambiguity; they then found a tendency for spin axes to 
be either parallel or perpendicular, with a shortage of pairs at /3 ~ 45°. 

All three methods outlined above for determining relative orientations have shortcom- 
ings. Determining the angle a between the spin vectors of two spiral galaxies requires high- 
resol ution imagin g or spatially resol ved s pectroscopy. In practice, t his has limited sample 



size; 



Heloul (119841 ). lOosterlool (119931 ). and IPestaiia fc Cabreral ( 120041 ) all had iV < 50 pairs 



in their samples. In addition, the assumptions must be made that dust lanes are on the 
outside of spiral arms and that spiral arms are trailing; these assumptions are not invariably 
true. Determining the angle /3 between the spin axes requires knowing only the apparent axis 
ratio q for each disk and the angle between their major axes; this can be found from lower 
resolution images, and permits larger sample sizes. However, because of the tilt ambiguity 
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of each disk image, the distribution of /3 is not, in the general case, determined uniquely. In 
addition, the assumptions must be made that disks are axisymmetric and that their intrinsic 
short-to-long axis ratio is known; these assumptions are not perfectly true. Determining 
the distribution of 0, the angle between the apparent major axes of binary spiral galaxies, 
can be done for even low-resolution images; smearing by the point spread function (PSF) 
will increase q for a galaxy, but will not strongly affect the position angle unless the PSF is 
severely asymmetric. This permits large sample sizes. However, even a perfect knowledge of 
the distribution does not yield a unique solution for /(a), the underlying distribution 
of angles between the spin vectors. 

The problem of determining the relative orientation of disks in binary spiral galaxies 
is a difficult one. Different investigators, even those using similar techniques or overlapping 
data sets, fiiid very different conclusions. These range from ra ndom orientation of spin axes 



(ISharp et al. 



parallel ( iHelou 



1979 



Poster 



1984 



Fhn 



oo 



19931 : ICervantes-Sodi et al.ll2010l ). to a tendenc y for axes to be 



19931 ) , to a tendency for axes to b e perpendicular (|So 



a tendency for axes to be either parallel or perpendicular (jPestana fc Cabrera 



ue 



1992), to 



20041 ). Our 



plunge into these troubled waters begins with selecting a sample of close pairs of disk galaxies 
from the Sloan Digital Sky Survey, as described in Section |21 Our initial analysis, described 
in Section [3] looks at the distribution of the angle between apparent major axes, for a 
relatively large (A^ = 747) sample of isolated pairs of galaxies. Then, in Section HI we add 
information about the apparent axis ratios q to find the distribution of /3, the angle between 
rotation axes. To eliminate the ambiguity in equation we look only at the N = 251 
isolated pairs for which at least one of the galaxies is nearly face-on [q > 0.9) or nearly 
edge-on (g < 0.3). Finally, in section |5l we discuss the implication of our results for the 
acquisition and evolution of angular momentum in binary spiral galaxies. 



2. DATA SELECTION 

We draw our sample of galaxies from the Sloan Digital Sky Survey Data Release 6 (SDSS 
DR6), which includes 9853 square degrees of photometric co verage and 7425 square degrees 



of spectroscopic coverage ( lAdelman-McCarthy et al.l l2008l ). The photometric imaging is 



done in 5 bands (ugriz); in this paper, we will use the imaging data in the r band. To 
build our sample, we start by selecting primary photometric objects that are classified as 
"galaxies" (obj_type = 3) in the r band and that have spectroscopic data with spectroscopic 
confidence parameter zConf > 0.35. We require that the redshift of each galaxy lie in the 
range 0.004 < z < 0.1; the lower limit on z eliminates contaminating foreground objects, and 
the upper limit reduces the possibility of weak lensing distortions of apparent shapes, and 
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eliminates in practice the necessity of applying K-corrections. To ensure that the galaxies 
in this sample are well resolved, we include only those images with r > 6.25rpsf, where r 
is the adaptive second-order moment of the galaxy image and Tpgf is the adaptive second- 
order moment of the PSF at the galaxy's location. The sample of galaxies satisfying these 
constraints contains = 166,853 individual galaxies. 

To extract a subsample of rotationally supported disk galaxies from our sample of 
166,853 galaxies, we start by using the SDSS parameter "fracDeV", which provides a mea- 
sure of the central concentration of a galaxy's light profil e. For each galaxy i mage, two 
models are fitted. One model has a de Vaucouleurs profile (Ide Vaucouleurslll948l ): 



I{R) = Jeexp (-7.67 [{R/Re 



a/4 



1]) 



(3) 



truncated beyond 7i?e to go smoothly to zero at 8i?e- The other model has an exponential 
profile: 

I{R) = hexp {-1.68 [R/ Re- 1]) , (4) 

truncated beyond 3i?e to go smoothly to zero at 4i?e. The SDSS data pipeline then takes 
the best-fitting de Vaucouleurs model and the best-fitting exponential model, and finds the 
linear combination of the two that best fits the galaxy image. The fraction of the total flux 
contributed by the de Vaucouleurs component is the parameter fracDeV, which is constrained 
to lie in the interval < fracDeV < 1. 

Previous studies using SDSS galaxies have used different cuts in fracDe V to separate 



l ate-ty pe ( low fracDeV) gal axies from early-type (high fracDeV) galaxies. IChang et al. 



(120061 ) and IShao et al.l (120071 ) , in their studies of late-type galaxies, found it useful to choose 
a sample with fracDeV < 0.5. In this paper, however, our aim is to produce a sample 
of strongly disk-dominated spiral galaxies. To t his end, we make a more s tringent cut in 
the fracDeV parameter, requiring fracDeV < 0.1. lUnterborn &: RydenI (120081 ) demonstrated 
that relatively luminous {Mr < —19) galaxies with fracDeV < 0.1 have the properties ex- 
pected of dusty, disk-dominated spiral galaxies: they lie primarily in the blue cloud of the 
color-magnitude diagram and they have colors and fluxes that are dependent on the appar- 
ent axis ratio q, with smaller values of q corresponding to redder, fainter galaxies. Moreover, 
visual inspection of SDSS "postage stamp" images reveals that g alaxies with Mr 5, —19 an d 
fracDeV < 0.1 appear to be primarily of Hubble type Sbc or Sc (jUnterborn fc RydenI |2008| ). 



The inclination-dependent dimming found by lUnterborn &: RydenI (120081 ) in their sam- 
ple of fracDeV < 0.1 galaxies is well fitted by the relation AMr = 1.27(logq')^. Because the 
dimming by dust is inclination dependent, a simple fiux-limited survey, such as the SDSS 
spectroscopic survey, will be biased against disks that are nearly edge-on (low q). To pro- 
vide a correction for this inclination-related bias, we require that all galaxies in our final 
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sample would be above the flux limit if they were seen edge-on. If an SDSS galaxy with 
fracDeV < 0.1 has q < 0.2, we assume that it is already edge-on. If it has q > 0.2, we 
compute its edge-on flux to be 

m,.(edge - on) = mr(obs) + 1.27[(log0.2)^ - (logg)^] . (5) 

To enter our corrected flux-limited sample, a galaxy must have mr(edge — on) < 17.77, 
corresponding to the completeness limit for the SDSS spectroscopic survey for galaxies. 

Finally, although disk-dominated spirals are we ll fitted by exponent i al pro files, so are 



other subclasses of galaxies, such as dwarf ellipticals. lUnterborn &: RydenI ( 120081 ) found that 
fracDeV < 0.1 galaxies dimmer than Mr ~ —19 have a m — r color that is independent of q. 
The lack of inclination-dependent reddening for the relatively low- luminosity galaxies reflects 
the fact that they are blue — r ~ 1.5) dwarf galaxies in which the stars and dust are not 
in orderly thin disks. To eli minate the dwarf galaxies from our sample, we start by following 



Unterborn fc Ryden ( 20081 ) in computing an approximate face-on absolute magnitude for 



each galaxy: 

M^(face - on) = Mr - 1.27(logg)2 . (6) 

To exclude the dwarfs, we require that Mr(face — on) < —19.4. Throughout this paper, 
when computing distances and absolute magnitudes, we assume a consensus cosmology with 
Ho = 70kms~^Mpc"\ Qm,o = 0.3 and Qa,o = 0.7. 

Applying the additional criteria that fracDeV < 0.1, (edge — on) < 17.77, and 
Mr(face — on) < —19.4, we reduce our sample size from = 166,853 galaxies to = 32,358 
luminous, disk-dominated spiral galaxies. To select binary spiral galaxies from this sample, 
we require that the redshift difference between two disk galaxies corresponds to a radial 
velocity difference Av < 300kms~^, and that the projected separation between them be 
Rp < I Mpc. To ensure that the binary system is relatively isolated, we also require that 
no other galaxy in the SDSS DR6 spectroscopic survey be within Av = 300kms~^ and 
Rp = 1 Mpc of either spiral galaxy in the binary system. These criteria leave us with 
Np = 747 isolated binary spiral galaxies. If we had chosen Av = 600kms~^ as our critical 
radial velocity difference (both for defining binaries and for defining isolated binaries), then 
we would have ended with Np = 699 isolated binary spiral systems, and the qualitative 
results of our study would have been unchanged. 

As a control, it will be useful to have a sample of luminous disk galaxies that are not 
part of a binary spiral system, and that are relatively isolated from bright galaxies of all 
types. To select a sample of isolated disks from our sample of A^ = 32,358 luminous, disk- 
dominated spiral galaxies, we require that no other galaxy in the SDSS DR6 spectroscopic 
survey be within Av = 300kms~^ and Rp = 1 Mpc of the spiral galaxy in question. These 



- 8 - 



criteria leave us with Ni = 16,814 isolated spiral galaxies. However, it is known that Sloan 
galaxies in close pairs have a si gnificantly different distribution in re dshift and stellar mass 
than non-paired Sloan galaxies ( lEUison et al.ll2008l : iPerez et al.ll2009l ). Thus, to ensure that 
the contro l galaxies are simila r in their properties to the galaxies in pairs, we adapt the 
method of iPatton et al.l (l201l[ ) for finding a matched control sample. We first go down the 
list of 2Np = 1494 disk galaxies in pairs, and find the isolated galaxy which best matches 
the redshift of each paired galaxy; once an isolated galaxy is matched to a paired galaxy, it 
is removed from further consideration. This gives a sample of 1494 isolated galaxies that are 
very closely matched in redshift to the galaxies in pairs. To increase the number of galaxies 
in our matched control sample, we iterate this procedure until the redshift distribution of 
the matched control sample deviates significantly from that of the galaxies in pairs. We 
found that 7 iterations gave a KS probability of Pks = 0.998 for comparison of the redshift 
distributions; going to an 8th iteration dropped the probability to Pks ~ 0.23. Thus, we 
adopt a matched control sample with Nc = 7 x 1494 = 10,458 galaxies. 

Figure [1] shows the cumulative distribution of redshifts for the disk galaxies in pairs 
(green line), for the matched control sample (red line), and for the complete unmatched 
sample of 16,814 isolated galaxies (black line). In this plot, the redshift distributions for 
the disks in pairs and the matched control sample are nearly indistinguishable (and in fact 
yield a KS probability of nearly one). Note that although we selected our matched control 
sample on the basi s of redshift, and not on both redshift and estimated stellar mass, as 



Patton et al.l ( 120111 ) did, we find that selecting purely on redshift gives similar (but not 



identical) distributions in color and absolute magnitude. 

The left panel of Figure |2] shows the cumulative distribution of g — r color for the 
disks in pairs (green line), for the matched control sample of isolated disks (red line), and 
the complete sample of isolated disks (black line). The control sample consists of galaxies 
that are redder than the disks in pairs, at a statistically significant level (Pks = 0.014); 
the median color of the control sample is g — r = 0.580, while the median color of the 
disks in pairs is g — r = 0.5 70. This difference in color may be due in part to starburst 
triggered in close encounters; iPatton et al.l ( 120111 ) found that "blue cloud" galaxies in close 
pairs {Rp < 80kpc, Av < 200kms~^) have g — r colors that are ~ 0.02 magnitudes bluer 
than a sample of galaxies, matched to have the same distribution of redshifts and stellar 
masses, that are not in close pairs. Since our sample of pairs contains primarily wider pairs 
(out to Rp = IMpc), it is unsurprising t hat o ur color difference, of ~ 0.01 magnitudes, is 
smaller than that found by iPatton et al.l ( 1201 ll ) . 

The right panel of Figure [2] shows the cumulative distribution of the corrected (face-on) 
absolute magnitude M/ for the disks in pairs (green line), the matched control sample of 
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isolated disks (red line), and the complete sample of isolated disks (black line). Although 
the complete sample of isolated galaxies is 0.148 magnitudes brighter, in the median, than 
the sample of paired galaxies, the matched control sample is 0.026 magnitudes fainter. For 
the remainder of this paper, we will compare the sample of Np = 747 disk galaxies in pairs 
to the matched control sample of isolated galaxies, and not to the complete sample of 16,814 
isolated galaxies, which skew toward higher luminosity, redder color, and higher redshift. 

The difference in median color between isolated disk galaxies and disk galaxies in pairs 
cannot be attributed to a difference in inclination-dependent reddening between the two 
populations. The green line in Figure [3] shows the cumulative distribution of the isophotal 
axis ratio q for the 1494 disk galaxies in pairs. The red line shows the distribution of q for 
the matched control sample of isolated galaxies. (If spiral galaxies were randomly oriented, 
infinitesimally thin, perfectly circular disks, then the distribution would be a straight diag- 
onal line from lower left to upper right on the plot; the lack of galaxies at small q indicates 
they are not infinitesimally thin and the lack of galaxies at large q indicates they are not 
perfectly circular.) A Kolmogorov-Smirnov test comparing the shapes of isolated spirals in 
the matched control sample (red) and the spirals in binaries (green) yields a probability score 
Pks = 0.367, consistent with their being drawn from the same parent distribution. 



3. ANALYSIS: APPARENT MAJOR AXIS ORIENTATION 

The first portion of our analysis involves studying the angle between the projected 
major axes of the two galaxy images in each of our Np = 747 isolated binary spiral systems. 
The angle is given by the relation 

COS0 = I cos(6'i — ^^2)1 5 (7) 

where 61 and 62 are the isophotal position angle for each galaxy in the pair, taken from the 
r band 25 mag arcsec"^ isophote. If the disks of the two galaxies are randomly oriented with 
respect to each o ther in three di mensions, then the distribution of (p will be uniform in the 



interval [0°, 90°] (IOosterlodll993l ) . The cumulative probability distribution for is shown as 
the black line in Figure HI the red line is the result for a distribution uniform in the interval 
[0°,90°]. A Kolmogorov-Smirnov test comparing the two distributions yields a probability 
score Pks = 0.892. Thus, the distribution of is consistent with the hypothesis of random 
disk orientations. 

A plot of the angle versus the projected separation Rp of the galaxies in each isolated 
binary is presented in Figure [51 The distribution of points reveals no trends with increasing 
Rp] in particular, close pairs have no tendency to be more strongly aligned. A KS test 
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comparing the distribution of (f) for the = 335 pairs with Rp < 0.5 Mpc to the distribution 
for the = 412 pairs with Rp > 0.5 Mpc yields Pks = 0.885, confirming the visual impression 
of Figure O 



4. ANALYSIS: ROTATION AXIS ORIENTATION 

Although a uniform distribution of is consistent with random orientations in three 
dimensions, it doesn't demand random orientations. In the general case, our knowledge of 
the distribution of (3, the angle between the spin axes of the galaxies, is impeded by the 
ambiguity of equation However, if we chose only those pairs which have at least one 
nearly face-on member, then we can assume an inclination cosii ~ 1 for one member, and 
thus 

/ 2_ 2X 1/2 

cos/3 COSi2 ~ ( Y , (8) 

where 7 = 0.22 is the intrinsic short-to-long axis ratio we use for the second galaxy. If the 
second galaxy has q2 < 0.22, it is assumed to be edge-on, resulting in cos/3 = 1. 

Similarly, if we chose only those pairs which have at least one edge-on member, then we 
can assume cosii ~ 0, and thus 

/ 1 - g2 \ 1/2 

cos /9 ^ sin ^2 cos ^ I 1 I cos(^i — 6*2) I . (9) 

VI - 7 / 

Identifying the spiral galaxies that are very nearly edge-on or face-on, using photometric in- 
formation alone, is necessarily approximate, given that spiral gal axies are not perfectly circu- 



lar and do not all have the same intrinsic short-to-long axis ratio. lUnterborn fc Ryderu (120081 ) 
fitted the distribution of intrinsic thickness 7 with a Gaussian distribution; for luminous 
fracDeV < 0.1 galaxies in the r band, they found a best fit of 7 = 0.21 6±0.067. Fitting the 



distrib ution of intrinsic disk ellipticity e with a lognormal distribution, lUnterborn fc Ryden 



(I2OO8I ) found \ne = —2.56 ± 0.91, implying a median ellipticity Emed = 0.077. 



For the purposes of this study, we assume that a galaxy with q > 0.9 is face-on, and a 
galaxy with q < 0.3 i s edge -on. If spiral galaxies truly had the shape distribution given by 



Unterborn fc RydenI (120081 ) and quoted above, and were viewed from random angles, then 
the galaxies falling into our "face-on" class would have (cos^ii) = 0.865 and the galaxies in 
our "edge-on" class would have (cos^ii) = 0.019. Of our 747 isolated spiral galaxy pairs, 
Nf = 94 pairs contain at least one face-on spiral and Ne = 171 pairs contain at least one 
edge-on spiral. Since 14 pairs consist of a face-on spiral mated with an edge-on spiral, the 
total number of pairs with at least one face-on or edge-on spiral is Np = 251. 
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In Figure [6l the black line shows the cumulative distribution function of cos /3 for the 
171 galaxy pairs containing at least one edge-on spiral; for these pairs, cos/3 is computed 
using equation The green line in Figure E] shows the cumulative distribution function of 
cos (3 for the 94 pairs containing at least one face-on spiral; for these pairs, cos (3 is computed 
using equation ([8]). Comparing the two distributions with a KS test, we find Pks = 0.823. 
The two distributions are consistent with their being drawn from the same parent population 
of cos (3, which is what we expect in the absence of inclination-dependent selection effects. 
If the spiral galaxies in binary systems were randomly oriented with respect to each other, 
then the distribution of cos /3 would be uniform in the interval [0, 1]. However, a KS test 
reveals that the distributions shown in Figure [6] are inconsistent, at a high level of statistical 
significance, with their being drawn from a uniform distribution of cos /3; for the pairs with 
an edge-on spiral, Pks = 0.021, and for pairs with a face-on spiral, Pks = 0.004. This 
is not necessarily an indication that the spiral galaxies are non-randomly oriented. The 
approximate values of cos/3 computed in equations ([8]) and (Q assume that spiral galaxies 
are perfect oblate spheroids, all with the same intrinsic axis ratio 7. These assumptions are 
known to be inexact. 

For a more useful test of the randomness of the orientation of disks in a pair, we compare 
the distribution of the computed values of cos (3 for the isolated pairs with the distribution 
of cosz for the = 10,458 isolated spiral galaxies in our matched control sample. The 
values of cos 2, computed using equation ([T]), contain the same inexact approximations as 
the values of cos/3, computed using equations ([H]) and In particular, the distribution 
of cos i for isolated spirals, given by the red line in Figure should be directly comparable 
to the distribution of cos /3 for pairs including a face-on spiral, given by the green line in 
Figure [6l in both these cases, we are computing the cosine of the inclination of a spiral 
galaxy relative to the line of sight to that galaxy. However, having a face-on spiral within 
Rp = 1 Mpc doesn't significantly affect a spiral galaxy's orientation relative to the line of 
sight; a KS test comparing the distribution of cost for galaxies in the control sample and 
coS;5 for pairs with a face-on galaxy yields Pks = 0.431. 

Our statistical tests are consistent with the hypothesis that the two disks in an isolated 
binary spiral system are randomly oriented with respect to each other. However, since 
structure on scales larger than ~ 1 Mpc has an influence on tidal torques, we also investigate 
whether the number density of galaxies on a larger scale may influence the distribution of 
cos (3 for binaries. For each of our 171 pairs with a face-on member and 94 pairs with an edge- 
on member, we determine the number of neighboring galaxies that are within a projected 
distance of Pp = 5 Mpc of either galaxy in the pair and with a radial velocity difference less 
than Av = 300kms~^ of either galaxy in the pair. To qualify as a neighboring galaxy, a 
galaxy must be in the Sloan spectroscopic sample, and its r band absolute magnitude cannot 
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be more than 1 magnitude dimmer than that of the lower-luminosity galaxy in the pair. The 
median number of neighboring galaxies, defined in this manner, was six. Consequently, we 
define pairs with six or fewer neighboring galaxies to be our "low-density" subsample, and 
pairs with more than six neighboring galaxies to be our "high- density" subsample. 

The left panel of Figure [7] shows the distribution of cos (3 for pairs in the low-density 
subsample, while the right panel of Figure [7] shows the distribution for the high-density 
subsample. In each panel, the distribution of cost for the isolated spiral galaxies in the 
matched control sample is repeated for reference (red line). If the simple picture of an 
isolated binary spiral system having parallel disks were correct, we might expect the pairs 
in the low-density subsample (left panel) to be show higher values of cos/3. In fact, we find 
that the distribution of cos /3 for both types of binary spiral systems (ones with edge-on disks 
and ones with face-on disks) is indistinguishable from the distribution of cosi for isolated 
spirals. The only marginally significant result we find, at the ~ 90% confidence level, is that 
in high-density regions (right panel of Figure [7]) binary spiral systems containing a face-on 
disk tend to have a distribution of cos /3 that is different from the distribution of cos i for 
isolated spirals: Pks = 0.088. The difference is such that the partners of face-on disks are 
more likely to be seen edge-on than an isolated disk would be. In fact, of the 42 partners 
of face-on disks located in high-density regions, 6 have q2 < 0.22. Of the isolated disks in 
the control sample, 4.9% have q < 0.22; given a probability P = 0.049 of an individual disk 
being edge-on, the expected probability of finding 6 or more edge-on galaxies in a sample 
of 42 would be P = 0.016, using the standard binomial formula. Although, in high-density 
regions, the partners of face-on galaxies have a high probability of being edge-on, the partners 
of edge-on galaxies do not have a comparably high probability of being face-on. 



5. DISCUSSION 

In summary, our study of disk galaxies in the Sloan Digital Sky Survey reveals no strong 
evidence for a preferred alignment of the individual disk galaxies in a binary galaxy system. 
For 747 close pairs of galaxies, the distribution of 0, the angle between the galaxies' projected 
major axes, is consistent with a uniform distribution between 0° and 90°. This is consistent 
with disks that are randomly oriented with respect to each other. 

Looking solely at 0, however, doesn't use all the available information. Computing 
cos/3, where /3 is the angle between the rotation axes of the disks, makes use of additional 
information - the apparent axis ratios of the two galaxies (equations [1] and [2]). Taking 
advantage of the loss of tilt ambiguity when one disk is face-on or edge-on, we can compute 
the distribution of cos/3 for these special pairs. Looking at 171 close pairs with an edge- 
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on galaxy and 94 close pairs containing a face-on galaxy reveals no statistically significant 
tendency for the two disk galaxies in a pair to be either parallel to each other (cos /3 ~ 1) or 
perpendicular to each other (cos/3 ~ 0). 

We find one curious result that is consistent with correlated orientations. In relatively 
high-density neighborhoods (more than 6 bright neighboring galaxies with Rp < 5Mpc), 
face-on galaxies within a close pair have an unusually high probability of having a edge-on 
partner; although just 4.9% of all disk galaxies in our sample have q < 0.22, 6 out of 42 
partners of face-on galaxies in high-density neighborhoods have q < 0.22. The physical 
significance of this statistically significant result is unclear. 

Obviously, the toy model for angular momentum acquisition, in which a pair of pro- 
togalaxies spin each other up and have antiparallel spin vectors and a tiny orbital angu- 
lar momentum, is inconsistent with our results. Acquisition of spin angular momentum is 
supplied by torques from an array of protogalaxies and protogroups in the vicinity. More- 
over, the initial angular momentum of a disk, acquired by tidal torques, can be modified 
by encounters and late infall. The continuing process of acquiring and modifying angular 
momentum appears, from our results, to produce binary spiral galaxies that have random 
relative orientations of their stellar disks. 
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Fig. 1. — The cumulative distribution of redshift z for the 1494 disk galaxies in pairs (green), 
for the matched control sample of 7 x 1494 = 10,458 isolated disk galaxies (red), and for the 
complete sample of 16,814 isolated disks (black). 
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Fig. 2. — Left panel: The cumulative distribution oi g — r color for the 1494 disk galaxies in 
isolated pairs (green), for the matched control sample of 10,458 isolated disk galaxies (red), 
and for the complete sample of 16,814 isolated disks (black). Right panel: The cuniulative 
distribution of M,., corrected to face-on using the presciption of lUnterborn fc Ryderu (j2008[ ) 
for the disk galaxies in isolated pairs (green), for the matched control sample of isolated 
disks (red), and for all isolated disks (black). 
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Fig. 3. — The cumulative distribution for the apparent axis ratios q of the 1494 disk galaxies 
in isolated pairs (green), for the matched control sample of 10,458 isolated disks (red), and 
for the complete sample of 16,814 isolated disks (black). 
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Fig. 4. — The cumulative distribution of the angle between the major axes of the galaxy 
images in the 747 isolated pairs of disk galaxies in our main sample (black); for comparison, 
a uniform distribution of is shown as the red line. 
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Fig. 5. — The angle between the major axes of the galaxies images in the 747 isolated 
pairs, plotted as a function of the projected separation Rp between the galaxy centers. 
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Fig. 6. — The green line shows the cumulative distribution of cos /3, estimated from equa- 
tion (jH]), for the Np = 94 isolated pairs in which at least one galaxy is face-on. The black 
line shows the distribution of cos/3, estimated from equation ([9]), for the A''^; = 171 pairs 
in which at least one galaxy is edge-on. A galaxy is judged to be edge-on if g < 0.3 and 
face-on if g > 0.9. For comparison, the red line shows the distribution of cosi, estimated 
from equation ([1]), for the matched control sample of isolated disks. 
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Fig. 7. — The cumulative distribution of cos /3 for pairs in the low-density subsample (left 
panel) and high-density subsample (right panel). Within each panel, the black line indicates 
pairs with at least one edge-on disk and the green line indicates pairs with at least one 
face-on disk. For comparison, the red line shows the distribution of cosz for the matched 
control sample of isolated disks. The values of cos /3 and cos % are estimated as described in 
the caption of Figure El 



